(DISCRETE) ALMANSI TYPE DECOMPOSITIONS: 
AN UMBRAL CALCULUS FRAMEWORK BASED ON osp(l|2) 
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Abstract. We introduce the umbral calculus formalism for hypcrcomplcx 
variables starting from the fact that the algebra of multivariate polynomials 
shall be described in terms of the generators of the Weyl-Heisenberg 
algebra. The extension of H\x} to the algebra of Clifford-valued polynomials V 
gives rise to an algebra of Clifford- valued operators whose canonical generators 
are isomorphic to the orthosymplectic Lie algebra osp(l|2). 

This extension provides an effective framework in continuity and discrete- 
ness that allow us to establish an alternative formulation of Almansi decom- 
position in Clifford analysis (cf. [381 1331 [36] ) that corresponds to a meaningful 
generalization of Fischer decomposition for the subspaces ker(D')'°. 

We will discuss afterwards how the symmetries of sl2(R) (even part of 
osp(l|2)) are ubiquitous on the recent approach of Render (cf. 1371 ). showing 
that they can be interpreted in terms of the method of separation of variables 
for the Hamiltonian operator in quantum mechanics. 



1. Introduction 

1.1. The Scope of Problems. In the last two decades considerable attention has 
been given to the study of polynomial sequences for hyperconiplex variables in dif- 
ferent contexts. For example, in the approach proposed by Faustino & Kahler 
(cf. [20]), rising and lowering factorials yield e.g. the classical Bernoulli and Eu- 
ler polynomials (cf. [IJ]) are the discrete analogue of homogeneous polynomials 
that appear in Fischer's decomposition involving difference Dirac operators. The 
hypercomplex generalization of this polynomials was studied recently in [35j by 
Malonek & TOMAZ in connection with Pascal matrices. 

Roughly speaking, the construction of hypercomplex Bernoulli polynomials shall 
be obtained via Appell sets O [34]. This was investigated at an early stage by 
Abul-ez & CONSTALES in [T] in terms of basic sets of hypercomplex polynomials. 
In this case, the resulting hypercomplex polynomials shall be interpreted as the 
Cauchy-Kovaleskaya extension of the rising factorials considered in [20] . For a fully 
explanation of Cauchy-Kovaleskaya extension we refer to [M] (Subsection II. 5); a 
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meaningful characterization of Cauchy-Kovaleskaya's extension in interplay with 
Segal-Bargmann spaces can be found in [10] (Subsection 2.2). 

Along with the construction of basic polynomial sequences in hypercomplcx vari- 
ables, other directions have been followed to construct Appell sets, namely by taking 
the Fueter-Sce extension of complex monomials z*^ (cf. |26j). using a Fourier series 
expansion of square integrable monogenic functions on the unit ball (cf. [5]) or 
alternatively using a Gelfand-Tsetlin basis approach (cf. [6]) that essentially is a 
combination between Fischer decomposition and Cauchy-Kovalevskaya extension. 

With respect to the discrete setting, it was further developed in the Ph.D thesis 
of Faustino (cf. [22]) that families of discrete polynomials can be constructed as 
a blending between continuous and discrete Clifford analysis giving an affirmative 
answer to the paper of Malonek & Falcao (cf. [34]). 

According to this proposal, discrete Clifford operators underlying the orthogonal 
group 0{n) were introduced by means of representations of the Lie superalgebra 
osp(l|2). Moreover, the refinement of discrete harmonic analysis follows from the 
representation of the Lie algebra 5(2 (R) as the even part of osp(l|2) while the 
blending between continuum and discrete Clifford analysis was obtained via a Shef- 
fer map (cf . [40] |39] ) that essentially maps the homogeneous polynomials onto basic 
polynomial sequences of binomial type. 

This approach combines the radial algebra based approach proposed by SOM- 
MEN in [43] with the umbral calculus approach postponed in [16] by Di BucCHi- 
ANiCO,LOEB & Rota. The main novelty of this new approach rests mostly from 
the fact that continuous and discrete Clifford analysis are described as realizations 
of the well-known Wigner quantum systems (cf. |45| ) on which the Sheffer map 
shall be interpreted as a gauge transformation that keeps invariant the symmetries 
of both systems (see e.g. [23] for a sketch of this approach). 

There were still to consider alternative constructions of discrete Clifford anal- 
ysis using different type of symmetries. One of them proposed in the preprint of 
Faustino & Kahler (cf. [21]), Clifford analysis on symmetric lattices corresponds 
to a mimetic description of Hermitian Clifford analysis on which the unitary group 
U{n) appears as the natural candidate to the induced representations for the algebra 
of Clifford- valued operators (cf. [7]). The major obstacle arising this construction 
follows from the fact that the multiplication operators XjT^'' and XjT^'' do not 
commute and hence there is no chance to get a radial algebra structure (cf. [43]). 
For a complete survey besides this drawback wc refer to [H] (Section 3). 

Quite recently, in the recent approach of De Ridder, De Schepper, Kahler 
& SOMMEN (cf. [T3]) the Wcyl-Hcisenberg symmetries encoded in the forward/ backward 
finite difference operators 9,^"' and multiplication operators x^T^"' were replaced 
by 'skew'-Weyl symmetries with the purpose to get, in analogy with the Her- 
mitian setting, linear independence between the vector multiplication operators 
= -^t ^^"^ ^ ^j=i ^1 '^J ' ^ result, the authors have shown 

that the Euler polynomials are the resulting discrete polynomials that yield a Fis- 
cher decomposition for + ~ ^^^h'' + ^J^h^ ■ Besides this approach 
there is an open question regarding the group of induced representations of such 
system. 

Let us turn now our attention for the Almansi decomposition state of art in 
Clifford and harmonic analysis. The theorem formulated below: 
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Almansi's Theorem (cf. [21 [3]) // / is polyharmonic of degree k in a starlike 
domain with center 0, then there exist uniquely defined harmonic functions on f2 
for-- ,fk-i such that 

fix) - foix) + \x\'f2ix) + ■■■ + 

corresponds to the Almansi decomposition for polyharmonic functions. 

One can find important appHcations and generahzations of this result for several 
complex variables in the monograph of Aronszajn, Creese &; Lipkin, [3], e.g. 
concerning functions holomorphic in the neighborhood of the origin in C". 

In the harmonic analysis setting, the importance of this result was recently 
explored by Render in [37], showing that for functions belonging to the real 
Bargmann space, there is an intriguing connection between the existence of a Fis- 
cher inner pair (cf. |37| ) with the problem of uniqueness for polyharmonic functions 
posed by Hayman in [27] (cf. [37] , Section 9) and as a characterization for the entire 
solutions of Dirichlet's problem (cf. [37], Section 10). 

In the Clifford analysis setting, the Almansi theorem shall be understood as 
a meaningful generalization of Fischer decomposition for hypercomplex variables 
without requiring a-priori a Fischer inner product (cf. [T3|, pp. 204-207). This 
result plays a central role in the study of polymonogenic functions likewise in the 
study of polyharmonic functions as refinements of polymonogenic functions. This 
was consider in the begining of 90's by Ryan to study the invariance of iter- 
ated Dirac operators in relation to Mobius transformations on manifolds. On the 
last decade Malonek & Ren established a general framework which describe the 
decomposition of iterated kernels for different function classes [33] |36] . Besides the 
approach of Cohen, Colonna, Gowrisankaran & Singman (cf. [12]) regarding 
polyharmonic functions on trees and the approaches on Fischer decomposition for 
difference Dirac operators proposed by Faustino & Kahler [20] and De Ridder, 
De Schepper, Kahler & Sommen [I3l, up to now there is no estabhshed frame- 
work on Almansi type decompositions as a general method for obtaining special 
representations for discrete hypercomplex functions. 

1.2. Motivation of this approach. The umbral calculus formalism proposed by 
Roman & Rota (cf. [321130]) have received in the last fifteen years the attention of 
mathematicians and physicists. Besides the papers of Di Bucchianico & Loeb 
(cf. [15]) and Di Bucchianico, Loeb, & Rota (cf. [H]) devoted to classical as- 
pects of umbral calculus, further applications were developed after the papers of 
Smirnov & TuRBiNER (cf. [42]) and Dimakis, Muller-Hoissen & Striker in 
the mid of the 90's (cf. [T^]) with special emphasis to the systematic discretization 
of Hamiltonian operators preserving Weyl-Heisenberg symmetries (cf. [32l [3T]). 
to the construction of Appell sets (cf. [4j) and complete orthogonal systems of 
polynomials (cf. ^7\) based on the theory of Sheffer sets likewise to the solution 
of the Boson-Normal ordering problem in quantum mechanics based on the inter- 
play between combinatorial identities using binomial sums with the construction of 
coherent states (cf. [3]). 

When we take the tensor product between the algebra of multivariate polyno- 
mials lR[a;] with the Clifford algebra of signature (0,n) in R", the resulting algebra 
of Clifford-valued polynomials is described in terms of Lie symmetries underlying 
the Lie algebra sl2(Il) and the Lie superalgebra osp(l|2) (see [HI [23l [11] and the 
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references given there) while the Fischer decomposition of the algebra of homo- 
geneous Clifford-valued polynomials in terms of spherical harmonics and spherical 
monogenics follows from the Howe dual pair technique (see [B] and references given 
there) applied to sl2(]R) x 0{n) and osp(l|2) x 0{n), respectively. The details of 
such technique can also be found in the papers of Howe (cf. [29]) and Cheng &; 
Zhang (cf.[S]). 

Although in the last years the Lie (super) algebra framework was successfully 
applied in Clifford analysis, these kinds of algebras are also ubiquitous e.g. in the 
old works of Wigner (cf. gS]) and TuRBiNER (cf. gTj): 

In [45] it was shown that the so-called Wigner quantum systems that describes 
a motion of a particle on the ambient space R" may be characterized in terms of 
symmetries of osp(l|2n); in [41] the eigenfunctions for the Hamiltonian operators 
were computed explicitly by taking into account the s[2(R) symmetries of such sys- 
tem while the eigenvalues were described as an infinite number of (unplaited) sheets 
lying on a Riemann surface. Quite recently, in |46| Zhang also apply this frame- 
work to the study of quantum analogues for the Kepler problem in the superspace 
setting. 

1.3. Organization of the paper. In this paper we will derive an umbral coun- 
terpart for the well known Almansi type decomposition for hypercomplcx variables 
by employing combinatorial and algebraic techniques regarding umbral calculus 
(cf. [mill]), radial algebras (cf. [43]) and the Howe dual pair technique confining 
nonharmonic analysis and quantum physics (cf. [28] [30]). 

We will start to introduce the umbral calculus framework in the algebra of 
Clifford-valued polynomials V := 'R[x] (E) Cio.n as well as the symmetries preserved 
under the action of the Shcffer map, showing that there is a mimetic transcription 
of classical Clifford analysis to discrete setting that generalizes complex analysis to 
higher dimensions (cf. [14[ I25j ). Roughly speaking, in umbral calculus the algebra 
of polynomials 'R[x] can be recognized as being isomorphic to the algebra gener- 
ated by position and momentum operators x'j and Oxj , respectively, satisfying the 
Weyl-Heisenberg relations 

(1) [O., , O.J = = [x;., .4], [O,^. , x',] ^ S.kid. 

Here and elsewhere [a, b] := ab — ba denotes the commuting bracket between a 
and b. 

Moreover, if we take ei, . . . , e„ as the Clifford algebra generators satisfying the 
anti-commuting relations {ej, 6^} := ej-efc -1-6^6^ = —26jk, umbral Clifford analysis 
(cf. [23] ) deals with the study of the algebra of differential operators 

A\g{x'j,Oj,.,ej : j = l,...,n} , 

For a complete survey besides this approach we refer to [22 (Section 3). 

Moreover, taking into account the action of the orthosymplectic Lie algebra of 
type osp(l|2) on the subspaces (2;')^kerD' we will derive some recursive relations 
(Lemma 13.11 and Proposition 13. ip and inversion formulae (Lemmata 13.31 and 13. 4p 
which allow us to decompose the subspace ker(D')'^ (the so-called umbral polymono- 
genic functions of degree k) as a direct sum of subspaces of the type Ps{x') kei D' , 
for s = 0,l,...,fc — 1, where Ps{x') stands a polynomial type operator of degree s 
satisfying the mapping property Ps{x') : kev{D'Y — > kerf. 
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This in turn gives an alternative interpretation for the results obtained by Ryan 
(cf. [38]), Malonek & Ren ([33l[36]) and Faustino & Kahler [20] in terms of 
the symmetries of osp(l|2). 

Finally, in Subsection 13 . 21 we will give an interpretation for the recent approach 
of Render (cf. [37]) showing that for a special choice of the potential opera- 
tor Vii{x') the Almansi decomposition encoded in the quantized Fischer inner pair 
{{2Vh{x'))'^ , (A')*"') is nothing else than a s[2(Il) based diagonalization of the Hamil- 
tonian H' = -^A' + Vn{x'). 

2. Umbral Clifford Analysis 

2.1. Umbral calculus revisited. In this section we will review some basic notions 
regarding umbral calculus. The proof of further results that we will omit can be 
found e.g. in [40l|39j|T5] or alternatively in [22] (see Chapter 1). 

In the following, we will set ]R[x] as the ring of polynomials over x — {xi ,X2,---,Xn) S 
R", by a = (ai, a2, ...««) the multi-index over Ng and by x" = x^^x'^'^ . . . x"" 
the multivariate monomial(s) over x. The partial derivative with respect to Xj will 
be denoted by ■= g|- while the n— tuple dx ■= {dxn dx2, ■ ■ ■ , dxr,) corresponds 
to the so-called gradient operator. 

Here and elsewhere, we will also consider the following notations: 

dl:=d^^d-^...dS:, al = aM...aJ., ^ ^ = |a|=E;=i«.- 

In addition we will denote by End(]R[a;]) the algebra of linear operators acting 
on ]R[x]. 

By means of the differentiation formulae dxX^ = for \a\ > \f3\ and dxxf' = 
(^p^'^y x^~" for \a\ < 1/3 1, it turns out the representation of the binomial formula 
in terms of dx' 

(2) fe+.)^=E (!!)^Y-"=E^H^^" 

|q|=0 ^ ^ |a|=0 

Linearity arguments shows that the extension of the above formula to IRfx] is 
given by f{x + y) = cxp{y ■ dx)f{x), where exp(y • dx) = E|^|=o Ir^x denotes the 
formal power series representation for the shift operator Tyf{x) = f{x + y). 

An operator Q S End(R[a;]) is shift-invariant if and only if it commutes with 
Ty = exp(y • d.x) for all P G R[x\ and y e R": 

[Q.Ty]P{x) Q{Ty_P[x)) - Ty_{Q P{x)) = 0. 

Under the shift-invariance condition for Q, the first expansion theorem (cf. |15] ) 
states that any linear operator Q : 'R[x\ ]R[ai] is shift-invariant if and only if Q 
is given in terms of the following formal power series expansion: 

oo 
\a\=0 

Set Ox = {Oxi,Ox2, ■ ■ ■ ,Ox„) as a multivariate operator. We say that Ox is 
shift-invariant if and only if Ox^, Ox2, ■ ■ ■ , Ox„ are shift-invariant too. Moreover, 
Ox is a multivariate delta operator if and only if there is a non-vanishing constant 
c such that Oxj {xk) = cSjk, holds for all j, fc = 1, 2, . . . , n. It can be shown that if 
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Ox is a multivariate delta operator, each Ox lowers the degree of P{x) 6 M.[x]. In 
particular O^^ (c) = for each non-vanishing constant c (cf. [35], Lemmata 1.1.8 
and 1.1.9) and hence, any multivariate delta operator Ox uniquely determines a 
polynomial sequence of binomial type {Va{x) : a G Nq}, (cf. [22] Theorems 1.1.12 
and 1.1.13): 

(3) ^M^ + y) = E ( f ) V^{x)V0.^{y) 

\P\=o ^ ^ 

such that Vq{x) = 1, Vq(0) = bafi and OxjVot{x) = ajVa-VJ{x)^ where Vj stands 
the j— element of the canonical basis of R". 

The Pincherle derivative of Ox with respect to Xj is defined formally as the 
commutator between Oxj and Xj-. 

0'xj{x) := [Ox,,x,]f{x) = Ox,{x,f{x))-x,{OxJ{x)). 

This canonical operator plays an important role in the construction of basic poly- 
nomial sequences of binomial type (cf. [THKITlll]). The subsequent results allows 
us to determine in which conditions (P'^.)~^ exists. 
We will start with the following lemma: 

Lemma 2.1. The Pincherle derivative of a shift-invariant operator Q is shift- 
invariant. 

Regardless the last lemma one looks to shift-invariant operators Q as formal 
power series (3(x) = ^if' obtained uiz the replacement of a; by i.e. t[Q(a;)] = 

Q{dx) where l : R[x] — >■ End(R[a;]) is defined as a mapping between the algebra 

of formal power series 'R[x\ and the algebra of linear operators acting on Il[a;]. 
According to the isomorphism theorem (see j39j . Theorem 2.1.1.), l is one-to-one 
and onto. This in turn leads to the following proposition: 

Proposition 2.1. A shift-invariant operator Q has its inverse if and only if Ql ^ 
0. 

From Lemma [2 . 1 [ notice that O'^^ is shift-invariant whenever Oxj is shift-invariant. 
Since from the definition O^^(l) = Oxj{xj) and Oxj{xj) is a non- vanishing con- 
stant, ProDOsition l2 . 1 [ asserts that (O'^. )~^ exists locally as a formal series expansion 
involving multi-index derivatives d". 

The former description in terms of Pincherle derivatives allows us to determine 
Va{x) as a polynomial sequence obtained from the action of (x')" := 11^=1 (^fe)"'"' 
with x'j. := Xfc(0^^)~^, on the constant polynomial $ = 1 (see also [39], page 51, 
Corollary 3.8.2): 

(4) v^{x) = {x:ri. 

The properties of basic polynomial sequences are naturally characterized within 
the extension of the mapping property ^'^ : i-^ Va{x) to R[x]. According to 
[39j . this mapping is the well-known Sheffer map that link two basic polynomial 
sequences of binomial type. It is clear from the construction that exists and 
it is given by the linear extension of ^'"^ : Va{x) H- x" to M[x]. In addition we get 
the following properties on R[x\: 

Ox, = 'i'xdx,'i'x\ and Xj = ^^Xj^'^^ 
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From the border view of quantum mechanics, the 2n + 1 operators a;^, . . . ,a;^, 
Ox^ , ■ ■ ■ , Ox„ and id generate the Bose algebra isomorphic to M,[x\ (cf. [IS])- Indeed 
for $ = 1 (the so-called vacuum vector) Oxj ($) = holds for each j = 1, . . . ,n 
while the raising and lowering operators, : Va{x) ^ Va+v^ (x) and Ox^ '■ Va(x) >-> 
C(j Va-vj{x) respectively, satisfying the Weyl-Heisenberg relations given by ([1]). 

Due to this correspondence, we would like to stress that the quantum mechanical 
description of umbral calculus give us many degrees of freedom to construct the 
raising operators Xj : Vq, (x) i— >■ Va+v^ (x) in such way that the commuting relations 
(U)) fulfil. In particular, in jl9i I17j it was pointed out the importance to consider 
the following symmetrized versions of Xj{0'^^)~^ 

(5) x'^ = lixAo'x,)-' + io'x,)-'^^) 

as a special type of canonical discretization. 

2.2. Basic operators. In what follows we will use the notation introduced in 
Section [2TTJ In addition we introduce C^o.n as the algebra determined by the set 
of vectors ei,e2,...,e„ satisfying the graded relations with respect to the anti- 
commuting bracket {a, b} = ab + ba: 

(6) {ej,ek} = ~2Sjk, 

The above algebra is commonly known in literature as the Clifford algebra of 
signature (0,ri) (cf. [Ij, Chapters & I; [25], Chapter 1) which corresponds to 
a particular example of an algebra of radial type. Indeed the anti-commutator 
{ej,efe} (scalar- valued quantity) commutes with all the basic vectors ej: 

(7) { [ej ,efc] ,e; } = 0, for ah j, fc, / = 1, . . . , n. 

For further details concerning the construction of C£o,n a-s an algebra of radial type 
we refer to [43] . 

Additionally, we will denote by "P = R[x} (8) C£o,n the algebra of Clifford-valued 
polynomials and by End('P) the algebra of linear operators acting on V. The 
Weyl-Heisenberg character of the operators x'j and Oxj combined with the radial 
character of the generators underlying Cio.n allows us to define umbral Clifford 
analysis as the study of the following algebra of differential operators: 

(8) Alg{a;;-,0,^.,e, : j = l,...,n} 

Furthermore, the umbral counterparts for the Dirac operator, vector variable and 
Euler operator, D' , x' and E' respectively, defined as follows 

(9) D' = Y,^,Ox,, 

71 

(10) x' = ^e,a;;., 

71 

(11) i?' = ^a;;0.,, 

correspond to linear combinations of the elements of the algebra defined in ([5]). 

In this context, the operators (P)- ([TT]) shall be understood as basic left endo- 
morphisms acting on the algebra End(7^). Along this paper we will use several 
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times the notation A' := X]j=i to refer the umbral counterpart of the Laplace 
operator A = Y.%i ^l^- 

The next lemma naturally follows from straightforward computations obtained 
by direct combination of relations ([T|) and 

Lemma 2.2 (cf. [12], Lemma 3.4.3, pp. 68). The operators x',D',E' G End{V) 
satisfy the following anti- commutation relations 

n 

{x',x'} = ~2"^{x'jf, {D',D'} = ~2A', {x',D'} = -2E' -nid. 

From the first relation of Lcmma[521 is a scalar- valued quantity while 

from the second relation A' := —{D'Y is a second order operator satisfying the van- 
ishing condition A'(l) = 0. On the other hand, in the third relation the action E' 
on End(7^) can be rewritten as a the following identity involving the anti-commuting 
relation between x' and D' : 

n 

(12) E' = J2x'fi^^ ^ --{{x',D'} + nid) . 

The action on both sides of the above identity on $ = 1 allows us to recast 
the dimension of the ambient space R" as n = —D'{x'l). So, the polynomial 
$ = 1 shall be interpreted as the corresponding ground level eigenstatc while the 
dimension of the ambient space R" appears as twice of the ground level energy 
associated to the harmonic oscillator containing n degrees of freedom. 

It is also clear from Lemma [2?2] that the anti-commutators {x' , x'}, {D', D'} and 
{x', D'} are scalar-valued. Thus, from (|9l),(IlO]) and ([12]) the operators {x' f, -{D'Y 
and E' shall be view as generalizations for the norm squared of a vector variable in 
the Euclidean space, the Laplacian operator and Euler operator, respectively. From 
the border view of quantum mechanics, the operators — ^(x')^ and —5 A' describe 
a spherical potential and the kinetic energy, respectively. 

Here we would like also to stress that the operator E' (see identity (ITT1) ) comprises 
at the same time the concept of directional derivative introduced by Howe (cf. |2H]) 
for quantum groups with the concept of non-shift- invariant mixed/number operator 
given by Di BUCCHIANICO, LOEB & Rota (cf. [16]). 

We will end this subsection by exploring and discussing some examples regarding 
the construction of the operators ([^- pip . 

Example 2.1. // we take Oxj — dxj, D' and x' coincide with the standard Dirac 
and coordinate variable operators, respectively: 

n 

while E' corresponds to the classical Euler operator E ~ X]j=i ^j^^xj ■ 

Furthermore, the continuum Hamiltonian i (—A + can we rewritten as 



-(-A+|xp)=-(Z?2-x2). 
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Example 2.2. Next we will consider D' as a difference Dirac operator given in 
terms of the forward differences d^^ f{x) ~ /(£+^vj)-/(a) g^ppg^-f^fj ^/jg /j^ . 

j=i i=i 

The square of D' corresponds to {D'Y — X]J=i(^2ftvj — "^Th-^- +id). On the 
other hand, the formal series expansion for d^^'' is given by d^'' = ^ (^cxp{hdxj) — id) , 
and [d^'' ,Xj] = T^vj = exp(/i9^^. ) corresponds to the Pincherle derivative for d^-' . 

Thus the operator x' corresponds to 

n n 

x' = E'^i XjT^hv, = x-jcxpi-hd^.). 
Alternatively, using relation the operator x' can also be taken as 

n 1 ^ 

x' = 2 E^'J +T-hvjXj) = 2^^^ {xj exp{-hda:^) + exp{-hdx^)xj) . 

Then we can consider two different constructions for E' (see relations ill]) and 

my- 

Hereby d^'' = (id — r_vj ) corresponds to the backward finite difference operator 
acting on the grid hZ. 

Example 2.3. Now we will replace the forward finite differences d^'' used to define 
D' in Example \2. 2\ by a central difference operator acting on hYl^ : 

0..f{x) = f(3i+h-^)-fi^-h^r,) 
The formal series expansion for these operators correspond to 

Ox J = ^ (exp {hdx^) - exp {-hd^.)) = ^ sinli {hd-x^) 
and moreover the formal series expansion for D' is given by 

1 " 

D' = ^E^j sinh(;ia,J . 

The square of {D'Y splits the star laplacian on a equidistant grid with mesh-width 
2h: 

Therefore, the construction of x' and E' shall be take into account the following 
formal series expansion for O'^. : 

O^./fe) = f(^+h.,)yi^^h.,) ^ ^^^^^^^^^^^ 
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Using the relation cosh(^hdxj) — cxp (— /i9a;^) (id — cxp (2/i9a;^)) combined 
with the standard Von Neumann series expansion of (id — exp (2hdxj ) ) ^ , we get 
the following asymptotic expansion for {0'^.)~^ : 

oo 

[O'xX^ = -h (id - exp {2hdx,))~^ exp {hd^^) = -/i^ exp {{2k + l)hdx^) , 

fc=0 

or equivalently {O'^. = —h X]fc°=o ^(2fc+i)/ivj • The above inverse only exists when- 
ever \\T2hv,\\ = \\exp (2hdxj)\\ < 1. 

Alternatively, we can express {0'^.)~^ using the following formal integral repre- 
sentation in terms of the Laplace transform (£/)(s) = /q°° e~''*/(t) dt (cf. [16j ); 

/* OO OO 

(O' )-' = -h e-^* exp{h{2t + 1)9,^ ) dt = -h e-''T^^2,+,^,^ dt. 
Jo Jo 

The umbral Dirac operator introduced in Exaniple l2.2l corresponds to the forward 
difference Dirac operator introduced by Faustino & Kahler in [20]. Here we 
would like to notice that in Example 12.21 the square {D')^ does not split the star 
Laplacian 



f{x + hv,) + f{x-hv,)-2f{x) 



which means that discrete harmonic analysis can not be refined in terms of discrete 
Dirac operators involving only forward differences (cf. PH]). 

As we see in Example 12.31 the computation of the inverse for O'^^ = cosh [hdxj ) 
is cumbersome and involves infinite sums or integral representations. However, in 
the case when periodic boundary conditions of the type x + hNvj = x for certain 
N G N are imposed on /iZ" (see [H], Section 5) it is possible to compute explicitly 
(0^.)~^ as a finite sum involving powers of T^vj ~ exp {hdxj)- 

On the other hand, contrary to Example 12.21 the operators x' and D' obtained 
viz the following intertwining properties on V: 

(13) ■^xD = D'^x. ^xx = x'-^x, ^xE = E'^x 

concern with the nearest neighbor points together with all the points contained in 
each direction hvj . Hereby, 'I'^^ is the Sheffer map introduced in Section 12.11 

Here we would also like to stress that = —{D')^ is supported on (2/i)Z". So, 
the periodicity as well as the coarsening of lattice is the price that we must pay in 
order to get discrete Clifford analysis as a refinement of discrete harmonic analysis 
underlying the orthogonal group 0{n). 

2.3. Orthosymplectic Lie Algebra Representation. The main objective of 
this subsection is to gather a fully description for the Clifford operators defined 
on Subsection 12.21 as a representation of the orthosymplectic Lie algebra osp(l|2). 
We will start to recall some basic definitions underlying the Lie algebra setting. A 
comprehensive survey of this topic can be found in [281 124j . 
The orthosymplectic Lie algebra of type osp(l|2) is defined as 

span {p^, p+, q} ©span{r^,r+} 
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equipped with tlic standard graded commutator [•,•] such that p ,p+,q, r and 
r+ satisfy the foUowing standard commutation relations (see e.g. [24]): 

[q,p±]=±p±, [p+,p-]=2q, 

[q,r±]=±ip±, [r+,r-] = iq, 

[p±,rT]=-r±, [r±,r±] = ±ip±. 

Here we would like to point out that on above construction, the Lie algebra s[2(Il) 
appears as a refinement of osp(l|2) in the sense that the canonical generators 
p~, p"*", q itself generate sl2(]R). In particular sl2(Il) corresponds to the even part 
of osp(l|2). 

In terms of the operators x', D' and E' -\- "^'id the subsequent lemma gives rise 
to an isomorphic characterization of osp(l|2) . We leave the proof of the following 
lemma to Appendix 1X1 

Lemma 2.3 (See Appendix [A|. The operators x' ,D' and E' + ^id generate a 
finite- dimensional Lie algebra in End['P). The remaining commutation relations 
are 

[x', [x'f] = 0, [x', -A'] = -2D', \e' + ^id, X-'] = x' 

[D\ [x'f] = -2x', [D\ -A'] = 0, [e' + ^id, D'] = -D' 

[[x'f, -A'] =a[e' + |id) , [E' + fid, -(x')^] = -2(x')^ [e' + ^id, -A'] = 2A' 

Furthermore, the standard commutation relations for osp(l|2) are obtained by 
considering the following normalization: 

p---iA', p+^-lixr-, q=i(^' + tid), r+ = ^zx', r- = ^zi?', 

and moreover, p"*" = p~ = and q = | {E' + ^id) correspond to the 

canonical generators of s[2(lR). 

In brief, the above description establishes a parallel with the continuum versions 
of Clifford analysis (cf. [14]) and harmonic analysis (cf. [30]) as representations of 
osp(l|2) and s[2(IR), respectively. This also establishes a link with the celebrated 
Wigner quantum systems introduced by Wigner in [JS] in the sense that the descrip- 
tion of the Clifford- valued operator in terms of the symmetries of osp(l|2) allows 
us to describe the motion of a particle confinining the quantum harmonic oscillator 
with n degrees of freedom. 

3. AlMANSI TYPE THEOREMS IN (DISCRETE) CLIFFORD ANALYSIS 

3.1. Main Result. In this section we will derive an Almansi type theorem based 
on replacements of the operators {x')'' by polynomial type operators Pk{x') such 
that the mapping {D')^Pk{x') : kerD' — ?> kerZ?' is an isomorphism. For a sake a 
simplicity, we leave for Appendix [B] the proofs of the technical results regarding the 
proof of the main result. 

We will start to pointing out the following definitions: 

Definition 3.1. Let Q be a domain in R" and k d JSl. A function f : 17 — > C^o.n 
is umbral polymonogenic of degree k if {L)')^ f{x) = for all x € f2. For k = 1 
f : 51 — > C^o.n is called umbral monogenic. 
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Definition 3.2. A domain fl C R" is starlike with center if for each x € 
txCzQ holds for any < t < 1 . 

For each k G IN, if the right inverse of [D')^ Pk{x') on the range of kerl?' exists 
we wiU define it by Q'^. : kerZ?' — > kerZ?', i.e. 

{D'fPk{x'){Q'J) = f, for all / € keri^'. 

Thus, the Almansi theorem can be formulated as follows: 

Theorem 3.1. Let be a starlike domain in R" with center 0. If f is a um- 
bral polymonogenic function of degree k in f2, then there exist unique functions 
/o, /i, . . . , /fe-2, /fc-i; each one umbral monogenic in such that 

(14) fix) = Poix')fo{x) + Pi{x')fi{x) + ■■■ + Pk-i{x)fu-i[x). 

Moreover, the umbral monogenic functions /o, /i, • ■ • , /fc-2, /fc-i ire given by the 
following formulas: 

(15) 

fk-i{x) = Q'k-iiD'f-'fisi) 
fk-2(M) = Qk-2{DT-'{id-Pk-i{x')Q',^,{D')''~')fU) 

fi{x) = QiD'(id- P2{x')Q2{D'f)---{id- Pi:^r(x')Qk^^{D')''-')f{x) 

fo{x) = (id - Pr{x')Q[D'){id - P2{x')Q'2{D'f) ■ • ■ (id - Pk-^{x')Qk-^{D'f~')f{x). 

Conversely the sum in |_?4[ ) with /o, /i, . . . , fk~2^ fk~i umbral monogenic in fl, 
defines a umbral polymonogenic function of degree k in Q. 



U' = 



Before proving Theorem 13. 11 we need some preliminary results. 

Lemma 3.1 (See Appendix [B|) . Let il be a starlike domain in R" with center 0. 
For any Clifford-valued function f(x) in fi, the following relations hold for each 
s G N; 

(16) D'iixyfix)) = ~2{xr-'u:j(x) + {-irix'YD'fix), 
where 

k id, if s ^ 2k 

-B' + (f + A:)id, if s = 2k + 1 ' 

From the above lemma, the next proposition naturally follows: 

Proposition 3.1 (See Appendix [B|) . The iterated umbral Dirac operator {D')^ has 
the mapping property 

[D'f : {x'YVcvD' {x'Y^'^keiD' 
for any s > k. Hereby, for each f{x) G kerZ?', 

(17) [D')' iix'Tfix)) = {-2Y{xr-'^UU^, . . . U:_,U'Jix), 
where the operators J7j are defined in Lemma \3.1\ 

Let be a star like domain with center 0. For any s > 0, we define the operator 
Is : G\n,C£o,n) C^{n,Cio,n) by 



(18) isfix)^ / f{txy-Ht. 

Jo 

In addition, we set Eg ~ sid + E. For s = we write E instead of Eq. 



1 

s-l. 
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Lemma 3.2 (cf. [33]). Let x G R" and ^l be a domain with ^ D [0,x\- U s > 
and f e C^(17, C£o,n); then 

(19) fix) = IsEJix) = EJJix). 

Sloppily speaking, the family of maps 1^ : C^(r2,C£o,n) — > C^{^7C£o,n) can be 
viewed as certain sort of right inverse for the operator Dx = fc=i ^j^kdxjXk in 
kerD. Indeed, if / is monogenic from Lemma [221 

D{xf{x)) = x{Dfix)) + D{xf{x)) = ~2Ef{x) - nf{x) = -2E^,^f{x) 

holds whenever Ox ~ dx- 

Finally, from Lemma [3?2l f{x) = -2E„/2 {-^In/2f{x)) = -\D{xI,^/2f{x)), 
showing that — ^/„/2 is a right inverse for Dx on the range kerZ?. 

On the other hand, when restricted to "P = ©fe^o^fc' where each G T'fe is a 
Clifford-valued homogeneous polynomial of degree k ( i.e. fkitx) — t''fk{x)), the 
family of mappings Is satisfy the equation 

1 

Isfkix)^ fkix)t''^-'-'dt= -—fkix), for all A en. 
Jo + s 

Hence Lemma [3.21 remains true for E' + sid = 'i'^^Eg'^'x in 'P, since it holds 
componentwise. 

Lemma 3.3 (See Appendix [B|) . There exists 1'^ : V ^ V such that 
(E' + sid)/^ = id = I'g(E' + sid). 
The next lemma will be also important on the sequel 
Lemma 3.4 (See Appendix [B|. If f eV, then 

(20) D'I'Jix) = I's+,D'f{x). 

For any k € Kq, denote by Q',^ = {-\f (f/fc)"HC^fe-i)"^ ■ • • (t^^)"^ where 

r i id, if s = 2k 

(21) {uT' = \ 

{ if.-2fc + l 

As direct consequence of \2Q\ . we find that I'gf{x) is umbral monogenic whenever 
f{x) is umbral monogenic. From the definition of Q'^. wc thus obtain 

(22) g;,(kcri:»') =keri:>'. 
Then the following lemma holds: 

Lemma 3.5 (See Appendix [BJ . For any umbral monogenic function f in Q,, 
{D'f [{x'fQ'Jix)] = /(£), x&^. 
Now we come to the proof of our main theorem for Pk{x') ~ [x')^: 
Proof of Theorem \3.1\ It is sufficient to show that 

haviD'f = kcr(D')''~^ + Pk-i{x')^CT D\ A: G N, 
where Pk-i{x') ~ {x'Y^^ . Notice that Lemma [3 . 5 1 states that 

(23) {D'fPk{x')Q't,^\A. 
We divide the proof into two parts: 
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(i) kerp')'' ^ ]icv{D'f-^ +Pk-i{x')keTD'. Since ker(L>')''"^ C ker(L>')'', we 
need only to show Pk-i{x')\ev D' C kcT{D')'' . For any g e kerD', by ([23]) 
and we have 

(I^O^Pfe-iCa;').?) = D'{{D'f-'Pu-^{x')Q',_M_,)-'g = D'{Q',_,)-'g = 0. 

(ii) kerp')'^ C kerp')''"^ + Pk-i{x')]<.CT: D' . 

For any / G ker(D')'^, we have the decomposition 

/ = (id - PkM^')Q'^-i{D'f-')f + Pk-i{x'){Q'k-,{Dr-'f)- 

We win show that the first sumniand above is in ker(L'')'^'^^ and the item 
in the braces of the second sumniand is in kerZ)'. Tliis can be verified 
directly. First, 

(D'f-\\d - Pk-i{x')Qk-i{D'f-')f = 
= {{D'f-^ - {{D'f-'Pu-i{x')Q'^_,){D'f-')f 
^ {{D'f-^ -{D'f-^)f ^Q. 

Next, since {D'f-^f<E kerL>' and Q;^,„i ker £>' C kerZ?', we have Q;^.„;^ V/ ^ 

ker D' , as desired. 

This proves that ker(D')'^ = ker(_D')'^^^ + Pk~i{x')kcY D' . By induction, we 
can easily deduce that ker(D')'' = kerD' + Pi(.T')kerD' + ^2(2;') ker Z?' + . . . + 
Pk-i{x')\iei D' . 

Next we prove that for any / e ker(L'')'"' the decomposition 

f = g + Pk-i{x')h, g e kcr(7^')'~\ fk e kciD' 

is unique. In fact, for such a decomposition, applying [D')^^^ on both sides we 
obtain 

[D'f-^f = {D'f-'g + {D'f-^Pu^^{x')h 

= iDr-'Pk^iix')Q'k^AQ'k^i)-'h 

Therefore fk = Q'k-i{D'f^^f^ so that 

g = f- Pk-i{x')fk = (id - P,_i(.t')QUiP')'"')/- 

Thus equations ([H]) and ([T5|) follows by induction. 

To prove the converse, we see from equation ^ of LcmmalXTjthat {D'Y+^{x'Ykev D' = 
{0} holds for any fc G N. 

Replacing k by j, we have 

{D'f{x'yVcvD' = {0} 

for any k > j. □ 
The proof of the above theorem can be interpreted as the following infinite 
triangle on which the subspaces ker(L'')'^ are despicted into columns. Each element 
of the triangle given by Proposition l3.1l corresponds to the action of osp(l|2) x 0{n) 
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on rows and columns: 



{0} 


kcrD' 


ker(L»')^ 


kcr(i:>')^ 


kcr(L»')^ 


{0} 




x' kcr D' ^ 




(a;')^kcr D' 




ID' 


ID' 


ID' 






{0} 


kcr D' ^ 


x' kcr D' ^ 


(x')^kerl?' 






ID' 




ID' 






{0} 


kerD' ^ 


x' kcr D' 


















{0} 


kerD' 



ID' 
{0} 



In those actions, the operator x' shifts aU the spaces in the same row to the right 
while the operator D' shifts all the spaces in the same column down. In particular, 
the {k + 1)— line of the above diagram corresponds to the action of (Z)')'^' on the 
subspaccs {x')^ kcr D' represented in (s + 1)— column. 

The next important step is the passage from the homogeneous operator of degree 
k, {x')'^, to a general polynomial type operator Pk{x') with the mapping property 
Pfc(x') : ker(D')'= ker D' . The corollary below gives a possible generalization for 
the construction of Pk{x'): 

Corollary 3.1. // Pk{x') = A'^. (x')'' + Rk{x') where A'^ is a Hilhert- Schmidt 
operator acting on V that satisfy the graded commuting relation [A'f.,D'] = akD' 
for some G IR and 

{D'f {Rk{x')f{x)) = 0, for all f e kcrD', 

then Theorem \3.1\ fulfils whenever the eigenvalues of the operator A'^, are greater 
than fcflfc. 

Proof. Starting from the definition of Pk{x') and using induction on k G No, the 
assumptions for A'j^ and Rk{x') lead to 

{D'f{Pk{x')f{x)) = {D'f{j^,{x'tf[x)) = {~kakid + ^^){D'f {{x'ff{x)) . 

whenever / belongs to kcrZ?' (i.e. / is umbral monogenic). 

From direct application of Proposition 13. 1[ the later equation becomes then 

{D')\Pk{x')f{x)) = {-2)\-kakid + AM ■ ■ ■ K-M{x)), 

where the operators U'^ are defined in Proposition 13. II 

Replacement of /(x) by S'J{x) = {-\f {U'M{U',_,)-^ . . . {U'^)-^f{x), on the 
above equation results in 

{D')\Pk{x')S'J{x)) = {-kakiA + j^^)f{x). 

Hereby (J/g)"^ are defined via equation ([2T|). 

Now it remains to show that — fcafcid + A'j. is invertible ensuring that Q'^. = 
S'^i-kokid + a;.)-! is a right inverse for {D')^Pk{x') : kerD' kerD'. 

If A'-^ is a multiple of id, Ok = and hence A'f. is invertible and the proof of 
CoroUarv 13.11 follows from Theorem 13. II 
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Otherwise, since A'j, is a Hilbert-Schniidt operator acting on V we conclude 
that A'^. has discrete spectra. Then, analogously to the proof of Lemma 13.31 (see 
Appendix |B]) A'j. is given by following the series expansion 

oo 

A'J{x) = Y,\k,sfs{x), 

where A^^s £ R correspond to the eigenvalues of A'j^. 

Thus — fca^id + A'^ is invertible whenever — fca^ + Aj. ^ is positive, that is A^.s > 

Finally, using the same order of ideas of the proof of Theorem 13. 1[ induction 
arguments lead to the following infinite triangle 
{0} kerD' 

{0} Po(a;')kerD' - 
{0} 



ker(D')^ 


ker(D')^ 


ker(L>')'* 


Pi(x')kerD' ^ 


P2(a;')kerL>' ^ 


P3(3;')ker D' 








Po(a;')kerD' ^ 


Pi(s:')kerD' ^ 


P2{x')kerD' 


\.D' 




ID' 


{0} 


Po(2;')kerD' ^ 


Pi{x')ker D' 






iD' 




{0} 


Po{x')ker D' 






ID' 






{0} 



This yields the following direct sum decomposition of ker(D')'^: 
ker(i:i')'' = ker(D')''"^ ©^A;-i(:E')ker£i' 

= ker(D')''^^^ ® ft-2(a;') kcrZ?' © Pfc-i(x') kerD' 

= Po (x') ker D' ® Pi {x') ker D' © . . . © Pfe_i (a;') ker D' . 

concluding in this way the proof of Corollarv l3.1l □ 
We will end this section by establishing a parallel between our approach and the 
approaches of Ryan (cf. [38|), Malonek & Ren (cf. [33l[36j) and Faustino &: 
Kahler (cf. HO]). 

Recall that Fischer decomposition ([Il|, Theorem 1.10.1) states the spaces of ho- 
mogeneous polynomials Vk are splitted in spherical monogenics pieces with degree 
not exceeding k: 

k 

^fc = II0^' CPk-snkcvD). 

s=0 

Moreover, from the mapping property given by Lemma l3.5l each Pk G Vk belongs 
to ker D^'^^ and hence from the intertwining property given by relations (|13p the 
Clifford- valued polynomial of degree k given by 'i>xPk{x) = Pk{x!)i belongs to 
ker(Z?')'^+^. Hence the following direct sum decomposition of ker(P'')'''"'"^: 

ker(P>')''^^ = ker D' © x' ker D' © (x')^ ker £>'©... © (x )'' ker £>', 

comprise the approaches of Ryan (cf. [SB]), Malonek & Ren (cf. [33]) (i.e. for 
'i'x = id) as well as the Fischer decomposition in terms of forward Dirac operators 
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obtained by Faustino & Kahler in [2^ if we consider the operators introduced 
in Example 1 2. 2 1 The fiexibihty of this approach allows us also to get the Fischer de- 
composition for several classes of finite difference operators like the finite difference 
operators considered in Example 12.31 

Remark 3.1. The replacement of [x')^ by the polynomial type operators Pk{x') 
given by Corollarv \3.1\ gives a parallel in continuum with the decomposition in terms 
of iterated kernels obtained by Ren &; Malonek (cf. |33j ] on which the operators 
Pk{x') shall be interpreted as quantizations of Clifford-valued polynomials of degree 



3.2. Parallelism with the Quantum Harmonic Oscillator. We will finish this 
section by turning out our attention to the quantum harmonic oscillator given 
by the following Hamiltonian written in terms of the potential operator Vh{x') = 



Hereby F' — —x'D' — E' corresponds to the umbral counterpart of the spherical 
Dirac operator (cf. [I4]). 

In order to analyze the 5(2(11) symmetries of "HJ^, we further introduce the fol- 
lowing auxiliar operator acting on the algebra of Clifford- valued polynomials V: 



The subsequent proposition gives a description of the Lie algebra symmetries un- 
derlying H'f^, showing that p+ = Vn{x'), = and q = f D' + + f id 
correspond to the canonical generators of sl2(Il). 
We start with the following lemma: 

Lemma 3.6. When acting on V , the operator F' commute with the operators E' 
and A'.' 



Proof. For the proof of [£",F'] = it remains to show that [E' ,x'D'] = since 
from definition [E\T'] = [E' , -x' D' - E'] = -[E',x'D']. 

From Lemma [Ql we get [E',x'] = x' and [E\D'] = -D' . This leads to 

E'{x'D') = {x' + x'E')D' = x'D' + {-x'D' + x'D'E') = {x'D')E', 

or equivalently [E',x'D'] ~ 0, as desired. 

In order to show that [F', A'] = 0, we recaU the relations [A', x'] = 2D', [A', E'] = 
2 A' and [A',D'] = that follow from Lemma [Ol This shows that 



k. 



^(^O'-l^' + f (r'-fid): 





[£:',F']=0, [A',F']=0. 



A' {x'D') = {2D' + x'A')D' 



2A' + {x'D')A', 



and hence, 



A'{x'D' + E') 



2A' + {x'D')A' + 2A' + E'A' = {x'D' + E')A' . 



Finally, taking into account the definition of F' the above equation is equivalent to 
[A', F'] = 0, as desired. □ 
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Lemma 3.7. When acting on V, the elements Vii{x') and J'^ are the canonical 
generators of the Lie algebra sl2(R). The remaining commutation relations are 



A' 



Proof. Recall that from Lemma 
canonical generators of sl2(R): 



P- = ¥ P 



A' 



A^ 
~2' 



a: - and q = are the 



[p , p^l = q, [q, p ] = p 



[q, p^l = q- 



Moreover [a;',p-] = -£»', = x' , [q, x'] = \x' and [13', p"] = = [a;',p+l 

Taking into account that 

Vn{x') = -p+ - fx' + ^ (F' - fid) and J^l = q + jD' , 



combination of the above relations with Lemma f3.6l results in the following identities 
in terms graded commuting relations: 



[-V',Vn{x')] 
[JLVnix')] 

This proves Lemma 13.71 



[P",P^JH 
q+|i3'; 



p-,_|.^' + i^(r'-|id) 



[q,p+]-f [q,x']-|[i3',p+]-(|)2[Z3',x'] 



= -P 



Vh{x'); 

[q> p"] 
p 



f (r' - fid) 



[f^',p- 



□ 



Proposition 3.2. The operators J^l^^T-L'j^ G End(V) are interrelated by the following 
intertwining property: 

n'f,exp{Vn{x'))exp (^-^^ = - exp (^^(x')) exp (^-^^ 

Proof. From Lemma [X71 the elements Vh{x') and correspond to the canon- 
ical generators of s[2(Il). 

From the above relations, it follows from induction over fc G IN that 

A' 1 



,Vn{x') 



leading to 
A' 



exp(F;,(.T')) 



kJL (yn{x')t-'' , [Ji Vn{x'f] = kVn{x'f-^ {Vn{x')f 



= j;,e^v{Vn{x')) , [Jie^v{Vn{x'))] = Vn{x')e^p{Vn{x')) 



Combining the above relations we get 



^ + Ji cxp {Vn{x'))\ = ( J-^ + Vn{x')) exp {Vn{x')) . 

is is equivalent to - Vn{x')^ exp {Vn,{x')) = exp {Vn{x')) + J"^ 
remains to show that + Jl^ exp (^4") ~ ^-'^P (~^) ^fi' 



This 
Not it 
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This statement is then immediate from the relation 



Y 



exp 



A^ 

T 



Therefore 

cxp (14(2:')) 

as desired. 



a; 



= -expiVnix')) 



A^ 

T 



Jn cxp - 



a; 

~2 



A' 

exp (Va (a;')) cxp ( — - 



□ 

We will finish this section by establishing a parallel with the recent approach of 
Render (cf. [37]). 

According to the definition of Fischer inner pair ([23: page 315) Corollary 13.11 
shows that when A'^, is a multiple of the identity operator, the pair {P2k{x'), (A')''') 
corresponds to a quantization of the Fischer inner pair that completely determines a 
(discrete) Almansi decomposition for polyharmonic functions (cf. [37] , Proposition 
20). Indeed, if for any umbral polyharmonic of function of degree fc on 17 (i.e. 
{A')''f{x) = holds on il) we take P2k{x') ~ {2Vh{x' j)'' it is straightforward from 
Corollary 13 . 1 1 that the following decompositions holds: 

fix) - Po{x')fo{x) + P2{x')fl{x) + ■■■+ P2k-2{x )fk-l{x) 

where /g, /i, . . . , /fe-2, /fe-i are umbral harmonic functions on 17 (i.e. A'/j(x) = 
holds for each j = 0, . . . , fc — 1 on ). 

Thus, it is also possible to obtain explicit formulae analogue to for umbral 
harmonic functions /o, /i, . . . , fk-i by considering the eyen powers of D' . 

Remark 3.2. Based on Lemma \3.7\ it is clear from the above construction that the 
functions /o, /i, . . . , fk-i obtained from Corollaru \3.1\ are solutions of the coupled 
system of equations: 

A'fk = 0, Jl^fk = (I + f )/fe. 

It is clear that in the limit ?i the above coupled system of equations approximate 
the umbral counterpart of spherical harmonics. In addition, from Proposition \ 3.S\ 
the composite action 0/ — exp (Vyi(a;')) cxp ^— on each fk span the eigenfunc- 
tions of T-L'f^. 

On the other hand straightforward computations combined Lemma \2.3\ leads to 
the following graded commuting property: 



E' 



-id 



,exp 



-D' 



-D'cxp -D' 



This yields the following intertwining property when restricted to the algebra V : 



E' + -id cxp 



-D' 



cxp 



D' Jk- 



showing that exp (^f maps the umbral harmonic polynomials of degree k onto 
umbral counterparts of spherical harmonics of degree k. 

In this case the action exp (—§£'') on V plays a similar role to the inversion of 
the Wick operator in Segal-Bargmann spaces underlying nilpotent Lie groups (cf. 

m)- 
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4. Concluding Remarks and Open Problems 

In this paper wc introduce an algebraic framework that can be seen as a com- 
prised model for Clifford analysis underlying the orthogonal group 0{n). This 
makes it possible to construct the associated operators and polynomials in the dis- 
crete setting starting from the equations and their solutions in continuum. The 
intertwining properties given by relations (fT3|) at the level of F,nd{V) gives us a 
meaningful interpretation of classical and discrete Clifford analysis as two quantal 
systems on which the Sheffer operator acts as a gauge transformation preserving 
the canonical relations between both systems. 

This approach can be interpreted as a merge between radial algebra approach 
proposed by Sommen [43] to define Clifford analysis with the quantum mechanical 
approach for umbral calculus described by DiMAKis, Hoissen & Striker (cf. [T9] ) 
and Levi, Tempesta & Winternitz (cf. [3T]). Based on the recent approach of 
Tempesta (cf . [33] ) we believe that this approach shall also be useful to construct 
polynomials in hypercomplex variables possessing the Appell set property. In this 
direction, the recent approaches of Malonek & Tomaz (cf. [35]) De Ridder,De 
ScHEPPER, Kahler & SoMMEN (cf.[l3]) and Bock, Gurlebeck, Lavicka & 
V. Soucek (cf. [6]) are beyond to the Sheffer set property. 

Here we would like to stress that contrary to the approaches of Malonek & 
Tomaz and Bock, Gurlebeck, Lavicka & Soucek on it is almost clear that the 
considered operators are generators of 5(2(11) (or alternatively s(2(C)) and osp(l|2) 
while the Appell sets are invariant under the action of the orthogonal group 0(n), 
in the approach of De Ridder,De Schepper, Kahler & Sommen it was not 
yet realized for which group the Appell sets (or more generally, the Sheffer sets) 
are invariant. 

Based on the recent paper of Brackx, De Schepper, Eelbode & Soucek 
(cf. [6]) and the preprint of Faustino & Kahler (cf. [21]), we conjecture the 
following: 

^All Hermitian operators represented in terms 0/5(2(11) 0,^6, osp(l|2) generators 
in continuum cannot be represented by 5\2{K) and 0Sp(l|2) generators in discrete 
but instead by quantum deformations of it '. 

For a nice motivation on this direction we refer to [21] (see Section 2) and also [H] 
(see Subsection 3.3) on which such gap was undertaken. 

In the proof of Almansi decomposition (Theorem 13.11 and Corollary 13. ip , the 
iterated (umbral) Dirac operators {D')^ play a central role. In comparison with 
[l4l [20] [22] [T3] we prove a similar result using the decomposition of the subspaces 
ker(D')'^ based on resolutions of osp(l|2) x 0{n) instead of considering a-priori a 
Fischer inner product. 

With this framework, Theorem 13.11 shows that the decomposition of ker(D')'^ in 
terms of osp(l|2) x 0{n) pieces yield the subspaces {x'Y kerD' for s = 0, 1, . . . , fc — 
1. Moreover, the replacement of {x')^ by a polynomial type operator Pk{x') in 
Corollary 13.11 gives an alternative interpretation for decomposing kernel approach 
proposed by Malonek & Ren (cf. [35]) as well as refines the Fischer inner pair 
technique used by Render in [37] to prove the Almansi decomposition in terms 
umbral polyharmonic functions. 

As it was observed along this paper the resulting approach based on represen- 
tation of the Lie algebra osp(l|2) as a refinement 5(2 (K) has a core of applications 
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in quantum mechanics that can further be consider to study special functions in 
Chfford analysis that belong to Segal-Bargmann spaces (see [11] and references 
therein). From the border view of physics, we have shown in Subsection 13.21 that 
the approach obtained by Render shall be described using a diagonalization in 
terms of s[2(Il). Indeed, Proposition 13.21 and Remark 13.21 explains the parallel be- 
tween the Almansi decomposition of the subspaces ker(A')'^ and the separation of 
variables of quantum harmonic oscillators (cf. [41 ] l42 l |46 ] ). 

One may further bring this technique in the future to construct new families 
of Appell/ShefFer sets for hypercomplex variables as well as to study Schrodinger 
equations on grids. At this stage, new families of discrete Clifford-valued polyno- 
mials like e.g. hypercomplex generalizations of Kravchuk polynomials (cf. [32] ) 
should appear. 
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Appendix A. Umbral Clifford Analysis 
A.l. Proof of Lemma 12.31 



Proof. Notice that the relations \D' , A' ] = = \x\ — (.t')^] are then fulfilled since 
{x'Y and —A' commute with all elements of End(7') (first and second relations of 
Lemma 12. 2[) . 

The proof of [E' + §id,.T'] = x' and [E' + fid,!?'] = -D' follow straightfor- 
ward from the Weyl-Heisenberg character of operators x'^ and Ox^ ■ Straightforward 
application of the above relations naturally leads to 



E'+'^\d,{x'f 



(x' + x' (e' + |id)) x' - x' (-x' + (e' + 



X -\- X 

2(x')^ 



n \ 
E' + -^d) 



[-D' + D' (e' + ^id)) D' - D' (^D' + E' + ^D'^ 
2A'. 
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Furthermore the relations [E' + fid, x'] = a;', [E' + |id, D'] = -~D' together 
with the third anti-commuting relation of Lemma 12.21 lead to 



D'{x'f - {x'fD' 



-2 
-2x' 



n . 
— ] 

2 



n . 
— ] 

2 



-2 ( + |id) - - D' (-2 (S' + |id 1 - a-'i?' 



-2 
21?' 



i?' + -id,Z?' 



-A'(x')' = i?' (-2.T' + (.t')'Z?') 

= -2W + (-2a;' + (x')'i?') i?' 

= -2{x' ,D'} - {x'fd.' 

= a(e' + |id) - {x'f^'. 



Finally, the combination of the relations [i?'+|id, x'] = x' and id, _D'] = -D' 

with the third anti-commuting relation of Lemma 12.21 leads to 



□ 



Appendix B. Almansi-type theorems in (discrete) Clifford analysis 
B.l. Proof of Lemma [3lTJ 

Proof. We use induction to prove (HH). Since {x', D'} = x'D'+D'x' = -2 {E' + fid) 
and Dg{x) = 0, we have 

(24) D'{x'g{x)) = -2 [e' + ^id) g{x). 

Next we show that, for any x G 51 and fc G N, 



^'((^')'^9(£)) 
i?'((x')''-'5(iE)) 



(25) 



~2k{x'Y^-^g{x)- 

-2(x')2('^-i) (£;' + (I + fc - l)id) g{x). 



This can be checked by induction. Assuming that (P5|) holds for fc. we shall 
now prove it also holds for fc + 1. We now apply the operator x' D' + D'x' = 
-2 [E' + fid) to the function {x'f^g{x): 

(26) x'D'iix'f'gix)) + i?'a;'((x')''5(i£)) = -2 (i?' + |id) ((x')'^g(ir)). 

By the hypothesis of induction, the first term in the left is equal to — 2fc(a;')^'°g(x), 
while the left side equals -2{x'Y^ [E' + {■ 



E' 



holds for all s > 0. 



s) id) x' -x' [E' 



2fc)id) g{x) due to the fact that 
slid) = E' +'^id,x' 



(DISCRETE) ALMANSI TYPE DECOMPOSITIONS 



23 



As a result, 

D'{{x')^f'+^g{x)) = -x'i?'(x'2'=5(x))-2(^' + fid)((.T')2'=5(x)) 

(27) = 2k{xy^''g{x) - 2{x')^'' {E' + (f + 2k) id) g{x) 

= -2(a;')2'= + (f + fc)id).g(s). 

This proves the second equality of The first equality of ([25]) can be proved 

similarly. This proves the identities ((25|) . □ 

B.2. Proof of Proposition [37TI 

Proof. In order to prove the above mapping property, we will derive (|17p using 
induction over fc G IN. First notice that for k = 1, relation ([TT]) automatically fulfils 
according to Lemma 13.11 

Next we assume that (fT7|) holds for any fc — 1, with fc e N. Hence, the action of 
D' on both sides of (fTT]) combined with Lemma [3.11 results in 

(^')'((^')V(3l)) = 

= (-2)'=+i(x')^-'=-iC/,:_,C/,:_fe+i . . . UUU'sfix) + {-ir-\xr-''D'gs-kix), 

with gs^kisi) = (-2)^C/,:_,.[/,:_,.+i . . . UUU'Jix). 

Now it remains to show that D'gs^k{x) — 0. If j is even, C/j = |id and hence 
[U'j,D'] = 0. Otherwise, from the second relations of [13] [£" + §,£)'] = D' com- 
bined with the definition of J7j for j odd results in [Uj, D'] = —D' . 

Thus, we have [U!j,D'] = — ^""'"^ ~^ D' for each j G IN and moreover for each 
g G kerD' the action of D' on Ujg{x) is equal to 

D'{U'^g{xj) = i— + U'^{D'g{x)) = 0, 

that is, D'U'jikerD') C kerZ?' for each j G IN. 

Finally, recursive application of the above relation leads to 

D'UUUU^, . . . U:_,U: (kerZ?') C kcr D' . 

and this results in D'gs-k{x) — 0, as desired. □ 

B.3. Proof of Lemma [Uni 

Proof. Take fkisi) G V such that E' fk{x) = kfk{x) holds for each k G INq. Hence 
for s > 0, the operator E' + sid has only positive eigenvalues of the form X = k + s 
which shows that the inverse of E' + sid tactically exists. 

For any f{x) = J2kLo fk{x) G 'P and s > 0, define : V ^ V &s the operator 
given by the series expansion 

oo 
fc=0 

Then we have 

oo oo 

{E' + s\d){I'J{x)) = E ((^' + = = /fe) 

fc=0 fc=0 
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Recalling the definition of in Section [2Jl the Clifford- valued polynomial yi>^ ^.fkj 
is homogeneous of degree k and hence 

leads to 

oo oo 

((£;' + Sid) fix)) = ^ /^((fc + s)Mx)) = ^ = fix). 

This shows that 1'^ is an inverse for the operator E' + sid, as desired. □ 
B.4. Proof of Lemma [131 

Proof. Starting from Lemma [231 we have —D' = [E' + ^id, D'], or equivalently, 

-D' = E'D' - D'E' = iE' + sid) D' - D' [E' + sid) 

by adding and subtracting (s — ^)D' on both sides of the first equation. This is 
equivalent to D' E'^ = E'^+iD' , where s ^ E'^ ^ E' + sid. 
Using the fact that E'^ ~ (J's) ^, wc end up with 

D'l', = I[,^,E',^,D'I[, = I's+iD'E'X - I's+iD'- 

□ 

B.5. Proof of Lemma [HH 

Proof. Denote gix) = Q'f^fix). From (|22p and the definition of Q'j., g is umbral 
monogenic in 17 and hence from equation (|17p of Proposition 13.11 we know that 

iDrHxTgix)) = i~2)''U[ . . . K_,U',gix). 

Thus iD'Y {ix')''Q'kfix)) ~ fix) follows directly from the above induced for- 
mulas. □ 
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